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Abstract 

The spin density wave instabilities in the quasi-one-dimensional metal (TMTSF)2C104 are stud- 
ied in the framework a matrix random phase approximation for intra-band and inter-band order 
parameters. Depending on the anion ordering potential V which measures the lattice doubling in 
the transverse direction, two different instabilities are possible. The SDWo state at low values of V 
is antiferromagnetic in b direction and has the critical temperature that decreases rapidly with V. 
The degenerated states SDW±, stable at higher values of V, are superpositions of two magnetic 
orders, each one on its subfamily of chains. As V increases the ratio between two components of 
SDW± tends to zero and the critical temperature increases asymptotically towards that of SDW 
instability for a system having perfect nesting and no anion order. At intermediate V, i.e. between 
the phases SDWo and SDW±, the metallic state can persist down to T = 0. 

PACS 75.10.Lp, 75.10.Fv, 74.70.Kn 



1 



Among the quasi-one-dimensional systems (TMTSF^X with the generic name Bech- 
gaard salts the material X=C104 has a particular status. If relaxed, it is a superconductor 
at ambient pressure and develops a field-induced spin-density-wave (FISDW) cascade for 
the magnetic fields B above 3 Tesla, well above the critical field for superconductivity. At 
even higher values of B it possesses a still controversial type of SDW ordering. If the 
sample is cooled more rapidly, the SDW order can persist even down to zero magnetic field. 
The critical temperature of this zero-field transition increases with the cooling rate up to the 
value of about 5.5 K for quenched samples. [0, || It is believed that this strong dependence 
of the SDW ordering on the cooling rate and its absence in carefully relaxed samples are 
related to the fact that (TMTSF) 2 C104 in its ground state acquires a long range order of 
anion molecules CIO4. This ordering sets up at the temperatures up to about 24K, inducing 
a doubling of the unit cell perpendicularly to chains, with the new lattice constant d = 2b, b 
being the inter-chain distance, as shown on Fig. |I| By varying the cooling rate one realizes 
a tunable unit cell doubling which in turn strongly influences the SDW ordering. 

The usual and the most direct measure of this doubling on the electron band properties 
is the band splitting V of a and (3 sub-families of chains (Fig. [l]).[f§] We take that V = 
— \V\ on a chains and V = \V\ on {3 chains. It was already pointed out that due to this 
splitting the two SDWs shown on Figs. [l](a) and |I|(b) are in competition. Jl], |5|, |6| Realizing 
that these SDWs cannot be treated independently (i.e. additively) as it was done in the 
existing bibliography, we develop a generalized matrix method within the random phase 
approximation (RPA). In this Letter we concentrate on (TMTSF) 2 C104 in the absence 
of magnetic field, and determine the phase diagram in whole range of values of the band 
splitting V, particularly for the physically relevant intermediate values, not properly covered 
in the previous analyses. 

As it is well known the Coulomb interaction between electrons leads to a SDW ordering 
via the mechanism of Fermi surface nesting. M If the chains were equivalent the nesting would 
lead to a SDW order with the magnetization modulation characterized by the wave-vector 
(2kp,7T/b) as shown on Fig. |l](a). One can imagine that in the limit of strong dimerization 
the nesting mechanism is again in action, but now separately for bands of a and (3 sub- 
families of chains from Fig. [I], and conclude that the critical temperature for each subfamily 
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FIG. 1: Magnetic pattern in direction perpendicular to chains for three characteristic cases of 
transverse SDW modulation: (a) the limit of weak anion potential V with p = n/b; (b) the limit 
of large V with p = ir/d; (c) V/tb moderate with p = ir/d transverse SDW order. 

is of the same order as that for equivalent chains provided the new nesting is again (close 
to) perfect. The new SDW then has doubled period 2c? = 46 in the b direction, as shown on 
Fig.|l](b) for the subfamily (3. Our analysis confirms that this indeed happens in the limit 
of large anion potential, V » tb, where t& is the interchain hopping integral. However for 
intermediate values of V the situation is more subtle due to the mixing of two SDWs. There, 
a new type of SDW dominates the response. It is a combination of two SDW orders on two 
different sublattices, each one forming an antiferromagnetic pattern in b direction, as shown 
schematically on Fig. [TJ(c) . We show that then the critical temperature is reduced or even 
completely suppressed. 

All modulations from Fig. [T] can be followed in terms of homogeneous and alternating 
SDW amplitudes denned on the dimerized lattice 

M h (x, n) = M a (x, nd) + Mp(x, nd + d/2) = -[M n (x, n) + M 22 (x, n)] 



M a (x,n) = M a (x,nd) — Mp(x, nd + d/2) = -[M 12 (i,n) + M 21 (x, n)}, (1) 

represented for later convenience in terms of SDW amplitudes on a{(3) chains M tt (^ = 
c +a(/3) ac - a il3)i an d i n terms of My = \f r + i <7\I/_ j (with i,j = 1,2) which are the SDW 
amplitudes in terms of the bond/antibond states 
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[cf, a (x, nd) ± c f:f3 (x, nd + d/2)], 



(2) 



with / = +,—. Here c± a ip\ are two component field operators for the right (left) Fermi 
sheet, and a is the Pauli operator in the spin space. Note that operators defined in the 



bond/antibond basis Mjj and NL h ^) depend only on the cell index n, while the local 
operators c and depend also on the position within the unit cell. In particular this 

means that c/ iQ (x, nd + d/2) = Cf^(x, nd) = and M a (x, nd + d/2) = M.p(x, nd) = 0. 

Below T c both and M a can attain non-zero mean values, with some specific common 
wave vector q = (k,p) and with the common orientation along some easy axis, say z. The 
staggered magnetization is then 

m z (x, R±) = (Ah ± A a ) cos + k)x + pnd] , (3) 

where the upper and lower sign stay for a (R± = nd) and j3 (R± = nd + d/2) chains re- 
spectively, and Ah/a are the mean values of the critical Fourier components (k,p) of M.h/ a - 
With A a ^ 0, Ah = and p = we get a simple antiferromagnetic order in b direction [Fig. 
|l](a)], like in the case of indistinguishable chains. We call this state SDWo. Another limit 
with doubled periodicity [Fig. |l|(b)] is realized for = A a or A^ = — A a with p = ir/d. 
In this case the SDW with antiferromagnetic modulation in the transverse direction can be 
realized either on the chains a or on the chains /3, each ordering having its own longitudinal 
modulation, with respective wave-numbers 2fcp + k\ and 2kp — k\. 

The intermediate ordering from Fig.|l|(c) is realized for \Ahj A a \ > 1. The corresponding 
ordered states with larger (2kp + k\) and smaller (2kp — k{) value of the longitudinal wave 
number will be denoted by SDW+ and SDW_ respectively. Fig. |](c) shows the ordering 
with Ah/A a < —1, i. e. with the larger SDW amplitude on (3 chains. Its longitudinal 
modulation has the wavenumber 2kp — k\. 

Written in terms of the Fourier transforms of the operators ([]]) in the perpendicular (b) 
direction with the Brillouin zone defined by (—ix/d < q < n/d), the interaction term of the 
Hamiltonian reads 

H i = ~tL [ dxY^[M[{x,q)M h {x,q) + Ml(x,q)M a (x,q)), (4) 

while the band contribution, written in terms of bond/antibond operators defined by eq.(^), 
reads 

TT/d 



p= — ir/d 



pd 

ivFp?,d x + 2r 3 t cos — + 2t' b cos pd — Vt\ 



*(x,P) ■ (5) 
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Here ^(x,p) = (^\(x,p),^\(x,p)), p's and r's are Pauli matrices in left-right and bond- 
antibond indices respectively, t' b is the parameter which introduces an imperfect nesting 
in the absence of the anion ordering, and 2N is the total number of chains. Also, being 
interested only in the SDW channel, we keep here only the forward part of the presumably 
local (intrachain) electron-electron coupling in Hj, omitting also the weak Umklapp term 
due to the existing slight lattice dimerization along the chain. 

Following the above qualitative arguments, it appears natural to construct the RPA 
formalism in the M aj /j basis. However, although the interaction term Hj is separated in M a 
and M/j, these SDW amplitudes are coupled through the non-diagonal part of the kinetic 
term ((5]), induced by anion ordering. Since the electron propagators in the bond-antibond 
basis are not diagonal, the RPA equations for M a ( ft ) amplitudes are also coupled. On gets 
the static SDW susceptibility in the form of a 2x2 matrix, 

(l-U Xaa)Xhh + U (Xha) 2 Xha 

Xha (l-U Xhh)Xaa + U (Xhaf 

(6) 

where indices i and j stand for h and a, Mi is the projection of Mj onto the spin quan- 
tization axis, D(q) = (1 - Uxhh)(l ~ Uxaa) - U 2 (xha) 2 , and Xhh,Xaa and Xah = Xha are 
corresponding bare (Hartree-Fock) susceptibilities. 
The diagonalization of the matrix leads to 



[x(q)] = ( (M ( (q)M?(q)>) = D~ l 



with 

1 

Xi/2 = 7^ 



(8) 



Xaa + Xhh ± yf (Xaa ~ Xhhf + ^(Xha) 

The SDW instability takes place for 

l-[/Xi(T c ,q c ) = 0, (9) 

q c being the wave vector at which xi (q) has the maximum. The corresponding ratio of two 
SDW order parameters from eq.@ is 



with 
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1 + 



Xa 



Xhh 



(11) 



\J (Xaa ~ Xhhf + 4:(Xha) 2 

The sign ambiguity in eq.fllQl) comes from the freedom of choice of the unit cell: (a, j3) or 

(P,a). 

It is convenient for further considerations to write the critical SDW susceptibility xi i n 
the form Xi(k,p) = [ln|? + Xres{k,p)], where Xres is always negative and contains all 
corrections to the pure logarithm. It can be derived after passing to the electron basis that 
diagonalizes the kinetic part (|5|) of the Hamiltonian. The result is 



+ 



1/2 



;i2) 



a 2 - P 2 )(a 2 p - (3 2 p )J2 XY ^xy) +lQUpJ2 X ^XY 

X,Y ' ^ X,Y 

Here (...) means the average along the transverse momentum (—ir/d < p' < tt/gT), \1/ is 
di-gamma function, X and Y take the values ±1, ^XY = + ^Qxy) w hh Qxy = 

v F k+2t' b [cos(p' + p) + cosp']-XA(p')-YA(p'+p) 



i 
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and A(p) = Jv 2 + [2i b cos 



pd 



a and (3 are the coefficients defining the unitary transformation that diagonalizes the band 
term (|), 

a(p') 



\ 



1 ± 



2t h \ cos(^ 



A(p>) 



(13) 



with the convention a p = a(p' + p) used to shorten writing in eq. (|12"1). The corresponding 
dispersions of two sub-bands are given by 



E ±f (k,p) = fv F k + 2t' b cos(pd) ± A(p) 



(14) 



The critical temperature is now the solution of the equation In ^% = Xres(T c , V, t b , t' b ; q c ) , 
where q c is always chosen to maximize Xres, and T c ° = T C (V — 0,t' b — 0) = eexp(— ^j^) is 
the tb-independent critical temperature for the asymptotic case of perfect nesting {t' b = 0) 
and vanishing anion order (V = 0); e is a constant of the order of the bandwidth. 

The wave vector dependence of Xres{k,p) is illustrated on Fig. ^]a, in which we chose 
t'b/tb — 0.1, V/tb = 1.2, and a very low temperature (T/t b = 0.0004) in order to make the 



whole fine structure of this dependence visible. There is one peak at qo = (4^/^,0), and 
two mirroring peak structures positioned approximately at q± = (/c ± ,7r/<i), with fc ± to be 
specified below. 

The peak at qo corresponds to the phase SDWo- It reproduces the standard SDW in- 
stability appearing in the absence of the anion ordering. [Tj Indeed it dominates over the 
peaks at q± for low values of V/fa, as it will be shown later in the discussion of the whole 
phase diagram. If in addition the temperature is not too much lower than t' b this peak is 
smeared and is simply shifted to k — 0, i.e. to the SDW with a pure 2k F longitudinal 
modulation. One gets the two-band version of the well-known imperfect nesting induced 
commensurate-incommensurate transition J7[, with rj w 1, i. e. with the SDWo order very 
close to an antiferromagnet in the transverse direction as shown on Fig. |I|(a). Here SDWo 
is stabilized as an interband process: the wave vector Q = 2k pd + q is the best nesting 
vector connecting the left Fermi sheet of the sub-band "+" with the right Fermi sheet of 



the sub-band whose dispersions are given by eq. (14). Note that in this range of wave 
vectors the anion potential V is also, beside the original imperfect nesting t' b , the source of 
deviation from the purely logarithmic susceptibility. It becomes clear from the expression 



(|i~2|) which after putting k = 0, p = and t' b = reads 



Xres(0,0) 



V 



2 



V 2 + [2t b cos(f 



1\ 1 i\V*+\2t b cos(^ 
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2 2 2ixT 



• (15) 



Note that one gets the same expression after putting k = 0, p = and t' b = into the 



result of Miyazaki et al ||11|| . For this particular case the off-diagonal element of the SDW 
susceptibility (^) vanishes, so that the matrix approach reduces to the standard scalar one. 

The peak at q + on Fig. 0a is in fact a plateau delimited on the axe p = ir/d by the wave 
numbers k± = (\/V 2 + Atf + V)/vp and k£ = 2^V 2 + 2t\jvF, as shown by the enlarged 
picture of this part of Brillouin zone on Fig.|^(b). This range of wave vectors is the intra- 
band nesting range for the sub-band E_ in eq. (|bj|) . Equivalently, the range q_ has the same 
role for the sub-band E + . The complete degeneracy of these two ranges would be raised 
after inclusion of higher corrections to the longitudinal linear dispersion in the one-electron 
spectrum assumed here. E. g. for the band-filling lower then one half the dominant SDW 
ordering is realized on the subfamily of chains which is higher in energy because it has a 



lower Fermi velocity and, consequently, a higher density of states. In that case, the minus 
sign is to be chosen in eq. (fL0|). 

The wave numbers kf and k^ delimit the region where after the translation by the wave 
vector Q + = 2kpa + q+ the left and right Fermi sheets cross each other, from the region 
where these sheets are not in contact. The difference Ak = kf — k£ measures the strength 
of effective imperfect nesting for the intraband SDW+ ordering. In the limit of large anion 
splitting, V ^> t b , it reduces to VF&-k rs t b /V 3 = 8e e g . In the language of the standard 
nesting model for SDW ordering J7| e e g is the effective next-nearest neighbor hopping in 
the diagonalized Hamiltonian, while t b takes over the role of an effective nearest neighbor 
hopping, i. e. it does not affect the perfect nesting for SDW+ and SDW_. This is in contrast 
to the interband SDWo ordering for which finite t b introduces an imperfect nesting, and tb 
keeps the nesting perfect. For temperatures T > e e g the fine plateau on Fig.|| is washed 
away and rounded to a single maximum at p = 7r/d and k\ ~ 2V(1 + tl/V 2 )/vF for V ^> tb- 

The magnetic pattern of the phase SDW_ is represented on Fig.[l|(c) for rj « 0.3. As 
V increases the parameter rj approaches the value 1/2 from below, i. e. the SDW_ mod- 
ulation in the transverse direction is closer and closer to the antiferromagnetic form from 
Fig.|I|(b). Furthermore, in the limit tb/V — > the nonlogarithmic term in eq.([T2|) reduces to 



its asymptotic value 



i. e. the critical temperature approaches 
:or the perfect nesting, T®. Note that in this limit Xres 
^(f) — ^^(1 + tt), as is seen from eq. (|I5D , i. e. on Fig. ||a the peaks at q± dominate over 
that at q . 

The phase diagram that results from the above considerations is shown in Fig.|3] in which 
the Coulomb coupling strength is parametrized by the perfect nesting critical temperature 
T c °. Fig.|3]a represents the case with t' b = 0. Then the low-V (SDWo) and high-^ (SDW±) 
phases are clearly separated by a "valley" in the range of intermediate values of V. The case 
with the imperfect nesting is illustrated on Fig.|3|(b) for t' b = 0.0375t b . 

The representatives of T C (V) dependences for the perfect nesting and highly imperfect 
nesting; CclSGS, clS well as the ^/-dependence of the parameter rj, are shown on Fig.^(c). The 
Coulomb interaction is chosen to give T c ° = 0.0108 x U b (i. e. T c ° = UK for t b = 300K). As 
was already pointed out, the critical temperature for SDW± phases, and the parameter r] as 
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FIG. 2: The quantity Xres 




FIG. 3: The phase diagram. Small discontinuity in the 77 vs V/4tb dependence is caused by the 
jump from k£ to (kf + fc^~)/2 in the Stoner criterion (0). 
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well [right side of Fig.^(c)] are not affected by t' b . On contrary, the critical temperature for 
phase SDWo is very sensitive on the parameter t' b , and is completely suppressed for t' b ~ T°. 

An attempt to make the correspondence of our phase diagram with the experimental one 
for (TMTSF) 2 ClO l i cannot be pursued without ambiguities. The main problem is a reliable 
choice of the value for the parameter t' b (beside the more or less well established values for 
other parameters, = 300.K" and T c ° = 13K. The value t' b = 11.25K is suggested from 
the experiments in quenched system M under the assumption that the sample is entirely 
free from anion ordering, i. e. that V = 0. Since one cannot exclude that even for the 
fastest quenchings some residual anion ordering with a finite, presumably small, value of 
anion potential V = Vres > remains, the above estimate gives an upper limit for t b in 
(TMTSF) 2 C10 4 . 

Choosing this value of t b , assuming that V = in the quenched samples, and expecting 
that V increases monotonously with the cooling rate (e. g. in measurements from Refs. 
0, ||), we come to the dashed line on Fig. |3|(c) as a fitting curve for (TMTSF) 2 C10 4 . This 
enables the estimation of the lower and upper limit for V in the relaxed salt. The latter follow 
from the width of the region where SDW phases do not exist, 0.1 < ^(relaxed ClO±)jt h < 
1.6. 

In conclusion, we have shown that the anion ordering in Q1D metals have some nontrivial 
consequences even in zero magnetic field. Fundamentally new concepts in this context are the 
need for a matrix SDW susceptibility in solving the interacting problem, and the subsequent 
introduction of the effective bare susceptibilities that diagonalize this matrix. Namely, all 
previous approaches to the interacting system with anion ordering used inadequately the 
Stoner criterion with simple SDW bare magnetic susceptibility || [HJ, or, equivalently, 
assumed that the order parameter is scalar in the bond-antibond space H |IT|| . Here we 
showed that the non-diagonal kinetic energy in bond-antibond representation is the one 
that dictates which effective bare susceptibility defines the Stoner criterion The novel 
proposition of the present work is the existence of the hybride SDW± phase at intermediate 
values of anion potential. Moreover, we have estimated precisely the lower and the upper 
limit of V for (TMTSF) 2 ClO A salt. 

The next important question is how the magnetic field influences the susceptibilities 
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Xi,2 given by eq. (§). Since the present matrix RPA is valid for finite magnetic fields as 
well, it remains to take properly into account the qualitatively more complex content of the 
kinematic part of the Hamiltonian (f|,|5]). We leave this part of work for future, expecting that 
the basic physics of the SDW phase in high magnetic fields is indeed that of two instabilities, 
as proposed already by McKernan et al 0). More specificaly, our analysis suggests that two 
relevant instabilities could be those appearing at wavevectors 2fci?a + q + and 2kpa + q_ on 
Fig. |. 
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